Tunneling in a uniform one-dimensional superfluid: emergence of a complex instanton 
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In a uniform ring-shaped one-dimensional superfluid, quantum fluctuations that unwind the order 
parameter need to transfer momentum to quasiparticles (phonons) . We present a detailed calculation 
of the leading exponential factor governing the rate of such phonon-assisted tunneling in a weakly- 
coupled Bose gas at a low temperature T. We also estimate the preexponent. We find that for 
small superfluid velocities the T-dependence of the rate is given mainly by exp(— c s P/2T), where 
P is the momentum transfer, and c s is the phonon speed. At low T, this represents a strong 
suppression of the rate, compared to the non-uniform case. As a part of our calculation, we identify a 
complex instanton, whose analytical continuation to suitable real-time segments is real and describes 
formation and decay of coherent quasiparticle states with nonzero total momenta. 

PACS numbers: 03.75.Kk, 03.75.Lm 
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I. INTRODUCTION 

Many one-dimensional (ID) systems share a universal 
low-energy description based on a complex order param- 
eter pj. Probably the most familiar example is a su- 
perfluid confined to a narrow channel, but — due to the 
well-known "duality" between bosons and fermions in one 
dimension — a similar description exists also for fermionic 
fluids. 

In accordance with the Bogoliubov-Hohcnberg theo- 
rem, no long-range order is possible in these ID systems, 
but the precise nature of fluctuations that prevent order- 
ing deserves a further discussion. At zero temperature 
(T = 0), perturbative fluctuations of the phase of the or- 
der parameter (phonons in a superfluid) cause a power- 
law decay of spatial correlations. At T ^ 0, the decay 
becomes exponential. However, a detailed study of the 
T = case reveals additional contributions to the corre- 
lation functions, of the form of a power law multiplied by 
an oscillating factor. 

While for weakly-coupled Fermi systems these oscillat- 
ing terms can be seen as a ID version of Friedel oscilla- 
tions, for weakly-interacting bosons their interpretation 
is not immediately obvious. It is possible, however, to 
interpret them as a consequence of nonperturbative fluc- 
tuations: instantons or quantum phase slips (QPS). The 
oscillatory dependence on the spatial coordinate can be 
traced to the fact (readily verified, see below) that each 
QPS changes the linear momentum of the superfluid com- 
ponent. 

The momentum production by QPS results from un- 
winding the order parameter and the corresponding 
change in the supercurrent. More formally, it can be 
viewed as a consequence of a special type of topological 
term, present in the action of a ID weakly-coupled Bose 
gas. We discuss this term in detail in the next section. 

A complementary picture is obtained by looking at the 
Lieb-Liniger spectrum (for a gas with a delta- function 
repulsion). Their results apply for periodic boundary 
conditions, i.e., ring geometry, which is the only case we 
consider here. In the limit L — > oo, where L is the length 



of the ring, the spectral branch associated with solitons 
H touches zero at momentum P = 2nn, where n is the 
gas density. This state corresponds precisely to the or- 
der parameter winding once as we go around the ring. 
Phonon-assisted transitions to this state will destroy su- 
perfluidity Our aim will be to calculate the rate of such 
transitions at low temperatures. 

Except for a brief summary in Sect. IIIII of results for 
QPS induced by a localized perturbation, we consider 
here only uniform ID superfluids. In the ring geometry, 
momentum in the longitudinal (x) direction is conserved. 
(More precisely, we should be talking about angular mo- 
mentum, but this distinction will not be important for 
our purposes.) Our goal was to see how momentum con- 
servation influences the QPS rate. Some results of this 
work have been presented in ref. Here we describe 
a different, more systematic method, which confirms the 
results of 0, but also allows us to obtain new results. 

At T = 0, the ID gas can be mapped on a two- 
dimensional (2D) model by introducing the imaginary 
(Euclidean) time r. QPS are vortices — or instantons — 
of this 2D model. Although this model is similar to the 
usual XY model, the topological term drives it into a dif- 
ferent universality class. The XY model, as the coupling 
is increased, undergoes a Berezinsky-Kosterlitz-Thouless 
(BKT) transition. In contrast, in the Bose gas, the cor- 
relation functions Q evolve continuously. We show in 
Sect. IIIII that, for a weakly-coupled uniform Bose gas 
at zero temperature, instantons and antiinstantons are 
bound in pairs by a linear, rather than logarithmic, po- 
tential. An extrapolation of this result to strong coupling 
implies that the breaking of instanton pairs, character- 
istic of a BKT transition, is not possible, a conclusion 
consistent with the expected Galilean invariance of the 
T = state. 

Instanton-antiinstanton pairs can unbind if there is an 
additional source of energy, besides the energy result- 
ing from unwinding the supercurrent. This possibility 
may be of interest for analog models of gravity 0,0. It 
has been suggested that one can use cold Bosegases to 
model cosmologically interesting spacetimes 0, U If one 
models an expanding spacetime by varying (decreasing) 



the coupling g between the atoms, as proposed in ref. 
8], then some of the energy supplied by this variation 
may become available for enhancement of QPS. We note 
also that a decrease in g causes the principal length scale 
of the gas — the "healing" length £ — to grow. Since the 
cross-sectional radius R of the channel is fixed, the ID 
regime £ > R can be reached. 

Another possibility, which is the main subject of this 
paper, is when the additional source of energy is a low 
but nonzero temperature. We discuss this case in detail 
in Sect. IIVI Thermally activated phase slips (TAPS) in 
the ring geometry have been considered, in the frame- 
work of the nucleation theory of Refs. 0, 0, [n|, m 
ref. [T^ . In the present paper, we concentrate on lower 
temperatures where, as we will see, the main mechanism 
for phase slips is thermally assisted quantum tunneling, 
rather than the over-barrier activation. One should note, 
however, that in the case of a uniform Bose gas, even at 
higher temperatures, the original LAMH theory 0, [n| 
probably needs to be modified, in order to account for 
momentum conservation. We discuss this point further 
in the conclusion. 

If the gas is uniform (or nearly uniform — the confining 
potential is smooth on the scale £, and its variations are 
small), QPS have a rather distinctive experimental sig- 
nature. Indeed, as we will see, unwinding momentum P 
from the supercurrent in such a uniform system is accom- 
panied by transferring a compensating momentum to the 
phonon "bath" . This is equivalent to creating a flow of 
excited atoms, which can in principle be detected experi- 
mentally. For example, we can start with a state with no 
supercurrent and let a QPS create a unit of supercurrent 
and a compensating normal flow in the opposite direc- 
tion (counterflow). The counterflow can be detected by 
the standard momentum imaging, i.e., opening the trap 
in one place. Note that no counterflow is expected for 
QPS induced by a highly localized perturbation, which 
breaks both momentum conservation and the Galilean 
invariancc (calculation of the rate for this case has been 
done in ref. 01). In that case, the final state of phonons 
has zero total momentum. 

However, perhaps the most immediate experimental 
consequence of the momentum balance during QPS in a 
uniform system is that it leads to a strong suppression 
of the rate, compared to the case of a localized pertur- 
bation The reason is that the momentum released 
by unwinding the supercurrent can only be absorbed by 
relatively high-energy phonons states, which are scarcely 
populated at a low temperature. 

From the theoretical perspective, one would like to 
understand in general how to compute the rates of in- 
stanton processes that transfer momentum between the 
background and excitations. The question is not limited 
to QPS in narrow superfluid channels but arises also in 
other contexts. For example, one can view the momen- 
tum transfer by QPS as a ID analog of the Magnus force 
in higher dimensions. This force acts on vortices moving 
in a superfluid and is known to suppress vortex tunneling 
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FIG. 1: A periodic configuration, in which instantons (open 
circles) and antiinstantons (shaded circles) are shifted relative 
to each other in space by Ax. At a finite temperature T, the 
period of the configuration is T = /3 = 1/T. 



|l4j . It is natural to ask if the suppression can be circum- 
vented by an inelastic process similar to the one we con- 
sider here. Furthermore, additional interesting physics 
emerges in cases when the order parameter is coupled 
to fermions, as in the case of BCS superconductors. The 
coupling to fermions opens a new channel of momentum 
production, due to the fermion zero modes at the instan- 
ton core. This channel is a ID analog of "momentoge- 
nesis" |15j by vortices in 2D arrays of Josephson junc- 
tions. Although in the case of ID superconductors the 
relevance of this channel is somewhat obscured by scat- 
tering of quasiparticlcs on the boundaries and disorder, it 
may still be of interest for interpretation of experiments 
jl6l Il7l| on superconducting nanowires. 

A much studied example of inelastic tunneling in field 
theory is instanton-induced scattering in gauge theories 
and their low-dimensional analogs . Indeed, the com- 
plex instanton that we will find in this paper is a gen- 
eralization of the (real) periodic instantons to the 
case when there is nontrivial momentum transfer from 
the background to quasiparticles. 

In ref. 0, we have identified periodic Euclidean con- 
figurations (of period (3 — 1/T, where T is the temper- 
ature) , which consist of chains of instantons and antiin- 
stantons shifted relative to each other by amount /3/2 
in the imaginary time and by some Ax in space, see 
Fig. n The tunneling rate has been obtained by in- 
tegratin g ov er Ax. Unlike the case without momentum 
transfer 19], this integration is nontrivial. Nevertheless, 
the leading exponential factor obtained in this way has 
a simple physical interpretation. It can be interpreted 
as the rate of tunneling between quasiparticle (phonon) 
states with momenta — P/2 and P/2, so that the change 
in the momentum of phonons precisely compensates the 
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momentum produced by unwinding the supercurrent. 

In the present paper, we derive this leading result for 
the tunneling exponent, and the first correction to it, in 
what we regard as a more systematic way. First, using 
a method developed in ref. [l9[ . we obtain the tunnel- 
ing exponent for a microcanonical state (fixed energy E). 
Then, we integrate the microcanonical rate over E with 
the Boltzmann factor exp(—E/T) to obtain the canonical 
(fixed T) rate. This method allows us to find directly the 
energies corresponding to the dominant phase-slip paths. 
In particular, we can show that at sufficiently low T 
thermally-assisted tunneling (as opposed to over-barrier 
activation) is indeed the dominant phase-slip mechanism. 

We also discuss in detail the instanton solution that 
saturates the rate at nonzero momentum transfer P. 
This instanton corresponds to a complex saddle point for 
Ax and is itself complex. Unlike the periodic instanton 
of ref. IS], it has no turning points. Nevertheless, it is 
possible to identify the initial and final states connected 
by this complex instanton and to reconstruct their real- 
time evolution. On the real-time segments, the solution 
is real and can be interpreted as formation and decay of 
coherent phonon states corresponding to the tunneling 
endpoints. A direct computation confirms that the total 
momenta of these states are —P/2 and P/2. 

II. THE TOPOLOGICAL TERM 

A weakly-coupled ID Bose gas can be described by 
the Gross-Pitaevsky (GP) Lagrangian, which in the Eu- 
clidean signature reads 

l e = tfw + ^-\d x ^\ 2 + % h/>| 4 - mI^I 2 • (i) 

2m 2 

Here r = it is the Euclidean time, m is the mass of 
the particles, g > is the coupling constant, and [i is the 
chemical potential. We assume that the system is subject 
to periodic boundary conditions in the x direction: tp{x+ 
L)=i>{x). 

Instantons are vortices of this theory in the (x, t) 
plane, corresponding to nontrivial winding of the phase 
of the order parameter ip. 

The number density of the gas is n — ip^ip an d can be 
written as a sum of the average density n and a fluctua- 
tion 5n: n(x,t) = n + 5n(x,t). We consider the uniform 
case, when n (the average) is independent of x, but will 
comment briefly on the effect of nonuniform n, such as 
resulting from a confining potential. 

At large wavelengths, fluctuations of the density are 
small, so ip has well-defined modulus and phase. We 
write tj) = in + tin) 1 / 2 exp(j#) and, expanding in small 
6n, obtain 

L B » i(n + Sn)d T 9 + ^-{d x 9f + Udnf . (2) 
2m 2 

Note that we impose no restrictions on the size of fluctu- 
ations of 0, in accordance with the absence of long-range 



order. Note also that in © we have omitted a term 
containing gradients of 6n. This is correct in the leading 
longwave approximation, but at shorter wavelengths that 
term will convert the purely acoustic dispersion law into 
the full dispersion of Bogoliubov's quasiparticles [2(J • 

The approximation J2J applies away from the instanton 
core, but not at the core itself. This is because the size 
of the core is determined by the healing length 

£ = (Agmn)- 1 / 2 , (3) 

while the long-wavelength approximation corresponds to 
wavenumbers k <C and so does not resolve the core. 

In this approximation, the density fluctuations are 
Gaussian, and can be integrated out in the path inte- 
gral. This amounts to using the equation of motion for 
5n and substituting the result, Sn — ~(i/g)d T 9, back into 
eq. (J2J). In this way, we obtain a phase-only theory with 
the Euclidean Lagrangian 

1 7) 

L E « md T 6 + ^(d T 6) 2 + i^{d x 6) 2 ■ (4) 

Instantons are singular solutions of this theory. The pres- 
ence of such singular solutions means that, even though 
the Lagrangian is quadratic in 9, the theory remains 
non-Gaussian. 

From eq. |Q} , we can read off the speed of Bogoliubov's 
phonons. We use units in which this speed is equal to 1: 

c s = (gn/m) 1/2 = 1 . (5) 

In these units, eq. © takes the form 

£ = Pffn)" 1 • (6) 

If it were not for the first term, the Lagrangian (0J 
would be that of the usual XY model. In that theory, in- 
dividual QPS can occur at T = in any current-carrying 
state, due to the energy released by unwinding the su- 
percurrent. This applies regardless of the strength of 
the coupling g. Moreover, at a sufficiently strong cou- 
pling (which is outside the domain of our semiclassical 
method) instanton-antiinstanton (IA) pairs unbind even 
in vacuum, resulting in a BKT phase transition. 

However, the first term in (@J drastically modifies the 
properties of the theory. We refer to this term as topolog- 
ical. It is of the same nature as the topological contribu- 
tion [2l| to the Magnus force acting on vortices in higher 
dimensionalities (2D and 3D). As we will now see, in ID, 
for Euclidean paths that contain instantons, the topolog- 
ical term gives a nonzero contribution to the Euclidean 
action. Being purely imaginary, that contribution can 
be interpreted as the interference phase between QPS at 
different locations. 

The simplest case when the the topological term can be 
seen to play a role is a single IA pair, with instanton lo- 
cated at Xq = (xo, To) and antiinstanton at x$ — (x' , Tq). 
Away from the cores, the density is approximately n, and 
the phase, in the absence of supercurrent, is 

0pair(a;, t) = arg[x - x Q + i(r - r )] - (x% <-> x%') . (7) 
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For periodic boundary conditions, this is an approxima- 
tion, which applies when all of the distances involved: 



\x^ — xfi\, \x^ — if j, and \xfi — xH \, are much smaller 



than the length L of the system. 

The time derivative of the configuration Q is 



X — Xq 



{x ~ x a ) 2 + (r - r ) 2 



(4 



(8) 



Integrating over r and x, we obtain the topological ac- 
tion: 



ASpan- = -2nin(x - x' ) 



(9) 



Since the variable conjugate to the position is momen- 
tum, we can interpret the action @ as resulting from 
production of momentum 



P = 2-im 



(10) 



by the antiinstanton and its absorption by the instanton. 

The above interpretation is of course nothing new and 
is readily confirmed by a direct calculation. Noether's 
momentum following from the real-time version of eq. 
Q is 



P 



dxtfr d x ip 



(11) 



The difference between the initial and final momenta can 
be written as an integral over a large closed contour in 
the (x, t) plane: 



Pr 



(12) 



where C runs clockwise. If the contour encloses a single 
instanton, without any phonon excitations, then away 
from the core the density is n, and Ijl2(l equals precisely 
— 27m. 



III. QPS RATE AT T = 

In general, the metastable states connected by the 
instanton are current-carrying, which implies that the 
phase 9 winds an integer number N times over the length 
L of the system. A suitable generalization of Q is then 



2tt 

3 (x,r) = 6» pair (a;,T) + -j-Nx 



(13) 



The second, winding, term corresponds to superfluid ve- 
locity 



V = 



2ttN 
mL 



(14) 



In what follows, we will always assume that V is much 
smaller than the speed of sound c s = 1 (although the 
more general case can be considered by similar meth- 
ods). For V <C 1, the field of the I A pair is only weakly 



affected by the superflow, and we can continue to use 
expression (JJJ. However, the presence superflow leads to 
an additional contribution to the action. To logarithmic 
accuracy, we obtain 



-iPAx + EAt - 



2tt d 
— lu- 
ff i 



O 



(15) 



where Aa; = x - xf , At = t - t^, d = (Ax 2 + At 2 ) 1 / 2 
is the instanton-antiinstanton separation, £ is the healing 
length |0, and 



E 



E N ^^^ = PV 
mL 



(16) 



Note that Em is precisely the energy released by a single 
instanton, as it unwinds the order parameter and reduces 
the supercurrent. 

The use of logarithmic accuracy means that we assume 
the separation d to be much larger than the instanton 
core size: d 3> £. This condition has to be verified a 
posteriori for typical configurations. 

To compute the rate of QPS at zero temperature, we 
need to integrate over all values of Aa: and At: 



1Z Im J dArdAx exp(—S pa i T ) 



(17) 



The rate computed in this way will be the inclusive rate, 
i.e., it will take into account the possibility of quasiparti- 
cle production in the final state. In other words, eq. l(T7|l 
can be viewed as an optical theorem for inelastic tun- 
neling. This relation between the imaginary part of the 
partition sum of IA pair and the rate of inelastic tunnel- 
ing [23| is familiar from the theory of instanton-induced 
cross-sections in particle physics [l^. Another way to 
look at it is to think about the IA pair as a bounce p3l | , de- 
scribing the decay of a metastable current-carrying state. 
Then, eq. (|17|l is the usual expression for the decay rate 
[24j], adapted to take into account the presence of the 
"soft" collective coordinates At and Aa; corresponding 
to the IA separation. 

As it is written, eq. I|17f) does not include the pre- 
cxponent and can be used to calculate only the leading 
exponential factor in the rate. (We will describe how to 
estimate the preexponent later.) Assigning the main role 
to the exponential factor implies the use of a semiclassical 
approximation and requires that the coupling g should 
be small. However, we will see that the main result — the 
absence of QPS in a uniform system at T = — can be 
plausibly extrapolated to larger couplings. 

The integral over Aa; is standard |25| : 



(a- 2 



cos(Pa;)da: 



P 



2IA7 



I>- 



-K v (P\At\) 



(18) 

where for our case 2v + 1 = 2n/g. Eq. (|18|) has to be 
integrated further over At with exp(— EAt). We see 
that this integral is convergent for E < P and divergent 
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for E > P. This means that at E < P the IA pair has a 
finite separation, that is, QPS are always bound in pairs 
and cannot occur individually. Note that the potential 
binding QPS at E < P is linear in |Ar|. 

Interpretation of the threshold at E = P is simple. 
Since instantons produce momentum, in a uniform sys- 
tem an isolated instanton has to be accompanied by 
production of quasiparticles (Bogoliubov's phonons) that 
carry that momentum away, so that momentum conser- 
vation is satisfied. But producing phonons with total 
momentum P requires, in a weakly-coupled gas, energy 
of at least E = P. Note that the condition of weak cou- 
pling is essential here. The coherent state of phonons 
that forms as a result of tunneling decays into individual 
quasiparticles, and in the weakly-coupled case we know 
their dispersion law — it is given by Bogoliubov's formula 
po|. So, we can explicitly verify that the energy of such 
a final state always exceeds its momentum (times c s ). 
This in general will not be true in a strongly interacting 
system (e.g., liquid helium). 

Now, if E is given by eq. (|16f) . i.e., the only source 
of energy is unwinding of the supercurrent, the condition 
E > P can never be satisfied in the superfluid state. This 
is obvious in the limit V <C 1, in which we have derived 
the pair action l)15|l . but the expression Ijltil) in fact holds 
also at larger V. We see that E > P implies V > 1, while 
according to Landau's criterion the superfluid state must 
have V < 1. Thus, in the absence of additional sources 
of energy, the superfluid state is always in the regime 
E < P, and at T = individual QPS cannot occur. 

Although we have obtained this result within the weak- 
coupling limit, the simple energetics underlying it allows 
us to speculate that it extends to arbitrary values of g. 
In other words, unlike the XY model, the uniform ID 
Bose gas has no BKT phase transition. This conclusion 
matches the observation that the correlation function of 
the Bose gas, found in ref. 0, evolve continuously with g. 
It is also consistent with the expected Galilean invariance 
of the T = state. 

The requirement for vortex unbinding, E > P, can 
be avoided in a nonuniform system, where momentum 
conservation need not be exact, for example, as a result of 
a short-scale perturbation. For a weak perturbation with 
length scale of order of or shorter than £, the exponential 
factor in the rate, to the same logarithmic accuracy, can 
be found by inserting a delta-function of Ax under the 
integral in i|17fl . so that 



2tt At 

7m, „ -Im / dArexpj -EAt In — 



(19) 



Eq. (|19|l coincides with the instanton rate that has 
been computed, for a microcanonical initial state of en- 
ergy E, in ref. [l^. That computation has been done 
in the context of the Abelian Higgs model, in the limit 
when it effectively reduces to the XY model. The salient 
point of the calculation is that the integral l|19|) is diver- 
gent (which, as we have seen, corresponds to unbinding 
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FIG. 2: Integration contour passing through the saddle point 

iH3J. 



of the IA pair) and has to be defined by analytical con- 
tinuation. The analytical continuation produces a finite 
imaginary part. (For a similar discussion in the context 
of dissipative quantum mechanics of a single degree of 
freedom, see Ref. 

Thus, we formally continue (|19fl to real-time separation 
At = -iAt: 



dAt exp -iEAt 



In- 



-At' 



(20) 



and then observe that the integrand has a saddle point 
at 



2tt 

At = At* EE I . 

gE 



(21) 



Deforming the integration contour so that it passes 
through the saddle point, as shown in Fig. [3 and replac- 
ing the exponent with the saddle-point value, produces 
the following exponential factor 19]: 



2vr 

K\ oc ~ exp \ — ln(gE(,) 
9 



(22) 



The exponent here has logarithmic accuracy, meaning 
that (|22|l applies only as long as gE£ <C 1. 

Although this section is devoted primarily to the T = 
case, we also list here for future reference the counterpart 
of eq. (|22|l for T ^ 0. This can be obtained from the 
action of the periodic instanton of Ref. by replacing 
the period with (3 — 1/T, or by integrating Ij22(l over E 
with the Boltzmann factor exp(— (3E). The result reads 



A? i ... ~exp<{ — ln(T£) 
9 



(23) 



and applies at T£ -C 1. In the context of cold Bose 
gases, perturbed by an external potential at a length scale 
shorter than £, this result was obtained by a different 
method in ref. 13], where in addition the preexponent 
was estimated. 

Because sharply localized perturbations that lead to 
these relatively large rates are unlikely to occur naturally 
in trapped atomic gases, it makes sense to inquire about 
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gentler sources of momentum non-conservation, such as 
a smooth variation of the trap potential. Since this ques- 
tion is somewhat outside the main subject of this paper, 
we will only address it qualitatively. Namely, we return 
to the action (|15fl . but now we do not assume that E and 
P in it are given by the specific expressions (|16(l and (|10ll . 
This describes qualitatively a situation when some of the 
momentum is absorbed by the potential. Note that for a 
smooth potential typical values of P will be close, even 
though not exactly equal, to the value 1)11) [1. Nevertheless, 
the regime E > P can now be realized. 

The same setup can be used to model the presence of 
an additional source of energy: we simply increase E rela- 
tive to ijlojt. One possible such source is time-dependence 
of the parameters, which, as discussed in the introduc- 
tion, is of interest for analog models of gravity. For ex- 
ample, decreasing the coupling g reduces the interaction 
energy, and one could imagine that some of the released 
energy becomes available for enhancement of QPS. Of 
course, there is no reason to think that the expression 
B15J1 with E > P will be literally applicable in this case, 
but one may hope that it will at least mimic some of the 
main features of the situation. 

As before, the integral over At is defined by analyt- 
ical continuation to real At = — iAr, so that eq. I|17(l 
becomes 



7T , Ax 2 - At 2 



U~ J dAtdAxe lPAx ~* EAt exp (~ In 



For E > P, the integrand has a saddle point at 



(24) 



/ jl a \ / » , , 2?r (E, P) 
(At, Ax) = (AU, Ax*) = i— J; ' . 25 

g E z — P A 

Deforming the integration contours so that they pass 
through the saddle point, cf. Fig. [21 we obtain the expo- 
nential factor in the rate as 



ft ~ exp - \n[(E 2 - P 2 )g 2 e 



(26) 



where the exponent again has logarithmic accuracy. 

We reiterate that, just as all the other rates computed 
in this section, eq. I)26|) is an inelastic rate — it corre- 
sponds to production of phonons with total momentum 
P in the final state. It can be viewed as a generalization 
of eq. (|22|l to the case when there is a nonzero transfer 
of momentum to quasiparticles. Eq. I|26|) has the ex- 
pected threshold at E — P, reflecting the requirement 
that production of phonons with momentum P takes 
at least E = P of energy (in units where the speed of 
phonons in equal to 1). At E < P, the integral in l|17fl is 
convergent, has no imaginary part, and the rate is zero. 
Note that the threshold is exponential: if the additional 
sources of energy and momentum are weak, the difference 
E — P is small, and the rate is strongly suppressed. 



IV. QPS RATE 

A. A summary 

We now turn to the main subject of this paper: calcu- 
lation of the QPS rate in a uniform system at a low, but 
finite, temperature T. The boundary problem, satisfied 
by the 5-matrix in the one-instanton sector at T =^ 0, is 
quite involved, and so we begin with a brief summary of 
the main points of the calculation. 

We consider only the case of low temperatures, 



T<^gn 



1 

2i 



(27) 



It would be interesting to extend the calculation to higher 
T . In the case of a sharply localized perturbation, we ex- 
pect a crossover to thermal activation at T ~ gn, i.e., 
when the period @ = 1/T of the periodic instanton be- 
comes comparable to the core size £, cf. ref. [l^. ^ n 
the Abelian Higgs model, such a crossover was indeed 
found numerically |27| . In a uniform Bose gas, where 
momentum transfer is necessary, it is not clear how such 
a crossover will occur and, in fact, even if it exists at all. 
In this case, the region T ~ gn can presumably be also 
addressed numerically, but in the present paper we will 
limit ourselves to a few comments on it in the conclusion. 

An intuitive approach to the problem at temperatures 
(I27|l would be to place the system on a cylinder of cir- 
cumference (3 = 1/T, i.e., consider configurations that 
are periodic in the Euclidean time r withperiod equal to 
f3. We have taken this approach in ref. [J. However, to 
make this approach rigorous, one has to prove that the 
configurations found in ref. [4| are indeed the dominant 
pathways for phase slips. A more systematic approach, 
which we take in this paper, is to first consider tunneling 
from a microcanonical state of a given energy E and then 
integrate over E with the Boltzmann weight exp(— (3E). 
This second method allows us to find directly the ener- 
gies corresponding to the dominant paths, and in partic- 
ular to show that in the regime l|27|) thermally-assisted 
tunneling is the dominant phase-slip mechanism, more 
important than over-barrier activation (if any). 

In addition, and curiously so, it turns out that to prop- 
erly explore the initial and final states connected by tun- 
neling, and also to compute the first correction to the 
semiclassical exponent, requires using a more precise dis- 
persion law for quasiparticles than the simple acoustic 
one. The method we use below allows us to take this 
modification into account without having to explicitly 
find the instanton solution corresponding to the more 
precise dispersion law. This method, developed in ref. 
[19j, is a sort of perturbation theory in energy — hence 
the limitation to the low-T range l)27|l . where the charac- 
teristic energies are not too large. In the present paper, 
we adapt this method to the case when instantons trans- 
fer momentum to quasiparticles. 

Now, the main difference from tunneling at T = 
is that at T ^ there are preexisting quasiparticles in 
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the initial state, and the tunneling path can make use of 
those. Indeed, we will show explicitly that tunneling now 
occurs between coherent quasiparticle states with total 
momenta —P/2 and P/2, so that the full momentum P 
is transferred but only energy E ps P/2 is required. This 
is in contrast to the T = case, where quasiparticles with 
total momentum P had to be produced, requiring energy 
of at least E = P. Accordingly, the leading exponential 
in the rate at T ^ is 

K T ~ exp(-P/2T + corrections) , (28) 

where the corrections are controlled by the small param- 
eter T/gn. Using the improved dispersion law for quasi- 
particles will allow us to obtain the first of these correc- 
tions. 



B. The boundary problem and the modified 
instanton 

The starting point point of our calculation is the ex- 
pression |l9j | for the microcanonical density matrix of 
phonons 

V E ~ 5(6 - E) (29) 
in the coherent-state representation: 

V E (b*,a)=Af J drjexp^-iEn + J dkb* k a k e %u,hl ^j . 

(30) 

The constant normalization factor Af is chosen so that 
Ve has the property of a projector: 

V E (b*,a)=V E (b*,a), (31) 

This makes it a projector on the subspace of a given 
energy E. 

The density matrix (|3L)fl will be our initial condition 
on a complex time contour shown in Fig. |3J cf. ref. 
|l9j |. The precise form of the contour will only become 
important later. For now, all that matters is that there 
is some initial Euclidean time Tj and some final r/, both 
of which can be complex but can be associated with the 
distant past and distant future, respectively. As typi- 
cal in problems where the evolution starts in the distant 

I 



Note that besides the usual path integral over the in- 
termediate values of the field 9 and the integral over rj, 
inherited from the projector H30JI . we have here also inte- 



past, the interaction is adiabatically switched off in the 
initial state, so that the Hamiltonian H in (|29|l is simply 
the Hamiltonian of free phonons. Accordingly, expiiukT]) 
factor in l|3U|l is a result of the free evolution of a coherent 
state | a): 

e i6r >\a k ) = \a k e 1 ^) , (32) 

where Uk is the frequency of the phonon mode with mo- 
mentum k. The coherent states are normalized by the 
condition (b\a) — exp(b*a). Combining this condition 
with eq. one can in fact quickly derive eq. (|3TJ|) . 





-r/4 




£ 

r/4 

X 


3 

Re? 





FIG. 3: Complex time contour for calculation of the S-matrix 
in the one-instanton sector. The circle denotes the instanton 
position. 



Note that we fix by hand the initial energy, but not 
the initial momentum. The most probable momentum 
of the initial state will be determined dynamically: it is 
the total momentum of the state J} fc \ak) corresponding 
to the most advantageous tunneling path. 

The probability of an individual QPS is 

Prob = Ti[SV E S*} = Ti[S E S E ] , (33) 

where S is the S- matrix, and S E = SV E ] we have used 
the projector property l|31|) . Tracing in eq. (|33|l is over 
all states that differ by unit of winding number from the 
supercurrent state on which the density matrix (|30|) is 
built. In other words, I133H is an inclusive probability. 

In the coherent-state representation, S E acquires a 
convenient form llfl: 



(34) 

I 

grals over the initial and final values of the field. These 
are needed to convert S E to the coherent-state represen- 
tation: expB, and expBf are the wave functions for the 



S E (b*,a) = J dd i de f V9dr]exp^-iEri + B i (a k e iw ^,e i )+B f (bt,ef)- J ' L E dr 
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coherent states \a k e tu ' kV ) and (b\: 



B t (a k e iu '"i,0i) = 
BM,6 f ) = 



dk 
dk 



-I afea _ fce 2 -<= - ^^(k^-k) + V^a^k^-^^ 



(35) 
(36) 



r 



where tildes denote spatial (in our case, one-dimensional) 
Fourier transforms, for example, 



dx 

—=i 
v/2tt 



(x,Ti)e 



ikx 



(37) 



The restriction that the winding number changes by 
one means that in 134(1 we consider paths of the form 

_ 2ir 

9(x,t) =d I (x-x a ,T-T ) + ^v{x,T) + —Nx , (38) 

where 9i is the one-instanton solution, and v is a fluctu- 
ation that has zero total winding (but may still include 
IA pairs). As before, we consider only relatively small N, 
when the superfluid velocity l|14|) satisfies the condition 
V<1. 

The approximate instanton solution, obtained from the 
longwave Lagrangian (@J, is 



r(x -x ,t- t ) = arg[x - x + i(r - r )] 



(39) 



We will see that at temperatures in the range 1)27(1 . the 
typical wavelengths of quasiparticles participating in the 
process are indeed large, much larger than the healing 
length £, so the longwave limit (Q is applicable. Never- 
theless, there are corrections to results obtained in this 
limit. In particular, for an accurate study of the saddle 
point that determines the QPS rate, the instanton l(39|) 
will have to be corrected ("modified"), to take into ac- 
count a deviation of the quasiparticle dispersion law from 
the purely acoustic one. 

Using 10} for Le in ((34(1 . we find that the integration 
over v{x, r) is Gaussian and gives rise to the free equation 
of motion: 



VV = 



(40) 



The integrals over the boundary values of v result in the 
boundary conditions (b.c.) 



Wfe5j(-fc) - Vi{—k) = 
u) k Vf(k) + Uf(k) = 



f%T ke ^(-T^v) ak 5 (41) 

^e^'bl . (42) 



Dots denote derivatives with respect to the Euclidean 
time r. Thus, the coherent-state parameters a and b* 
determine, through the b.c. ((41|) and 142(1 . the non- 
Feynman parts of the fluctuation v — the negative fre- 
quency part at Tj and the positive frequency part at r/. 



The Fourier transform of the instanton field l(39|) itself 
is computed at real To, such that Rer^ < tq < Rer/, 
and real xq and then analytically continued to arbitrary 
complex values. We find 

— 6i(k,n) = -^ =e UJ "^- T ^ +lkxa R{k) , (43) 

Vs v 2 ^fc 



-u>k{Tj—To) — ikxo 



where 



1/2 



R{k) , (44) 



(45) 



Because the instanton field satisfies Feynman b.c. in 
both directions, it does not directly participate in the 
b.c. (flTTl . (14*2*1) . However, as seen from (|3l*|) . ([3*H|l . it acts 
as a source for the coherent-state parameters a and b* . 
In the expression ((4511 . 



u k = \k\ , 



(46) 



in accordance with the fact that the solution ((39(1 was ob- 
tained from the longwave limit |(3}, in which the phonon 
dispersion is a simple acoustic one. As it turns out, im- 
portant phonon momenta in our case are of the order 



fc^(^ T) -l/3 ln l/3| 



(47) 



where T 3> £ is the period of the configuration (see be- 
low). Since this momentum is much smaller than £ , the 
acoustic approximation 1(46(1 is adequate for use in the co- 
efficients R(k). However, as will be shown below, in the 
exponentials exp(±cjfcr) we need to use a more precise 
approximation to Bogoliubov's full dispersion law |20|: 



u; k = (k 2 +k"e) i/2 ^\k\+^k\ 3 e 



(48) 



This means that the instanton that we will be using is in 
fact modified relative to the simple expression ((39(1 . This 
modified instanton could in principle be obtained from 
the Fourier transforms ((43(1 , 144(1 , in which we substitute 
the more accurate dispersion law l(48|) in the exponents. 
One of the advantages of the present approach, however, 
is that it will allow us to obtain many interesting physical 
quantities without ever needing the explicit form of the 
modified instanton. 
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Microcanonical rate 



Gaussian. After some algebra, we obtain 



The solution to the boundary problem H40l) - (|42|) is 



v(x,t) 




a k e 



— T-\-irj)-\-ikw _i_ l* uj^t— ikx 



(49) 



where v' is a solution to l|40[) with Feynman b.c. A non- 
trivial v' is only possible because the region of applica- 
bility of H4Ufl has a "hole" at the instanton core. To find 
z/, we need, in principle, to consider the equation for the 
entire fluctuation of the order parameter, including both 
the modulus and the phase. In other words, instead of 
(13811 we would write 



tp = (i/ji + Sip) ex-p(2iriNx/L) 



(50) 



where ipi is the instanton solution, and Sip is a fluctua- 
tion. Instead of the longwave limit (QJ, we would have 
to consider the full Lagrangian . The integration over 
Sip will no longer be automatically Gaussian, but it will 
become such in the leading semiclassical approximation. 
The corresponding equation for the fluctuation is 



V 2 Sip = , 



(51) 



where T> 2 is the second-order differential operator in the 
instanton background. Fortunately, in what follows we 
will not need the explicit form of (|51|l . but only the gen- 
eral properties of its solutions. 

There are two types of solutions to eq. ifST]) . Solutions 
satisfying the Feynman b.c. are the zero modes of the 
operator T> 2 , associated with the collective coordinates 
xq , To of the instanton 1|39[1 . The coefficients with which 
these zero modes occur in Sip are arbitrary and need to 
be integrated over. These integrations can be converted 
in the usual way into integration over the collective co- 
ordinates. 

The other type of solutions are delocalized modes. 
These contain non-Feynman parts, so the coefficients 
with which they appear in Sip are fixed by the b.c. The 
main idea of the perturbative method developed in ref. 
|19| is that when the typical momenta k of excitations 
involved in tunneling are small, the delocalized modes 
can be approximated by the plane-wave solution given 
by eq. (|^9l with v' = 0. We stress that this "pertur- 
bative" method is not an expansion in small coupling g. 
Instead, corrections to the plane wave, which are due to 
scattering of the plane wave on the instanton core, are 
controlled by the parameter According to the esti- 
mate (|47f) . this parameter is small. Thus, to the leading 
order, we can simply neglect v' and use in (|34[) the field 
(|5"5jl with v given by the plane wave only. 

Next, we observe that the tracing in the expression 
(|33|l for the probability can be done in the coherent- 
state representation by integrating over a and b with 
exp(— a* a — b*b) as the measure, and this integration is 



Prob- 



where 



dxodtodx' dt' di]dr/' e 



-2S +iPAx-iE N At-iEQ+W 



(52) 



W = J dk\R\ 2 Z(k)< 



,iukC\ e -iKAX , JKAX-iuj k <: 



2]- 

We have introduced the following notation: £ = 77 — 77', 
Ax = xq — x' , At = to — t' 0l KAX = ujkAt — kAx, and 



E(fe) = (1 - e^)- 1 . 



(54) 



R is given by eq. H45|) . As in eq. (|24|l , we have regularized 
the divergent integrals over To and Tq by continuation to 
real time. 

The doubling of the number of integrations in (|52() has 
to do with the presence of two path integrals in (|33|) : one 
in Se, and the other in S E . If Se is associated with an 
instanton, then S E can be associated with an antiinstan- 
ton. The first three terms in the exponent of (|52|l are 
the sum of the actions of a single instanton and a single 
antiinstanton. In particular, to logarithmic accuracy, 



S = - In — , 

9 6 



(55) 



where the size of the system L is used as an infrared 
cutoff. Comparing (|55|) to the action (|15fl of an instanton- 
antiinstanton (IA) pair, we see that L now appears in 
place of the IA separation d. This is because we are now 
computing not the action of a pair, but the sum of the 
individual actions, each of which is infrared-divergent. 
As we will see, the dependence on L will be removed 
by the W term (|53ll . whose presence reflects the non- 
vacuum nature of the initial and final states of phonons. 
This is similar to how the IA interaction at T = reflects 
a non-vacuum final state, a correspondence well-known 
from the studies of instanton-induced cross-sections in 
particle theory [H llH . 

Nontrivial integrals in (|52|l are those over the relative 
positions At, and Ax. The remaining integrals, those 
over the "center-of-mass" positions, simply produce pow- 
ers of space and time volumes, which are either absorbed 
in the normalization of the projector Ve-, or factored out 
when we compute the probability per unit length and 
unit time, i.e., the rate. 

We begin by integrating over £ and Arj. These integrals 
can be done by steepest descent, and the correspond- 
ing saddle-point conditions have simple physical meaning 
In our case, the saddle-point conditions are 

E = 2 J dkv k \R\ 2 e iUkC (cos K AX - 1)S 2 , (56) 
E N = J dkuj k \R\ 2 {e lKAX _ e <"*C-«fAX) S _ ( 57 ) 
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These can be rewritten as 



dkcu k a k a k 



E N = dkbJk{b* k b k - a* k au) 



(58) 
(59) 



where a, a*, &, and b* are the saddle-point values cor- 
responding to the Gaussian integration that led to 1)52(1. 
These values characterize the most probable initial and 
final states and will be discussed in detail in the next 
subsection. For now, we will only need the correspond- 
ing quasiparticle densities: 



a* k a k = \R 



b* k b k = \R 



2 sin 2 [w fc At/2 - fcAcc/2] 

sin 2 KC/2) ' 
2 sin 2 [w fe (C - At)/ 2 + fcAcc/2] 
sin 2 (w fc C/2) 



(60) 



(61) 



where the upper sign corresponds to a, and the lower 
one to b. Notice that l(58|l is a natural expression for 
the energy of the initial state, while 159|) expresses en- 
ergy conservation: the change in energy of the phonon 
subsystem equals the energy produced by unwinding the 
current. 

It is of interest to consider both the case when Ace is 
real (which is its original domain), and the case when 
Acc is purely imaginary (which is where the saddle-point 
for it will be found). In either of these cases, the saddle 
points for £ and At are purely imaginary, in particular, 



C = iT 



(62) 



with T > 0. As we will see, T is the period of the 
configuration. 

In this paper, we consider only the limit when the 
energy En released by unwinding the current is much 
smaller than the typical thermal energy E (although the 
more general case can be considered similarly). Then, 
the left-hand side of (|57|> can be set to zero, and we find 
that 



At= C -= l -T 
2 2 



(63) 



These are the same saddle-point relations as those found 
in ref. 19], but obtained here for a somewhat more gen- 
eral situation — when an instanton causes a nonzero mo- 
mentum transfer. 

The energy condition (|56|) can now be rewritten in the 
form 



E 



1 



dkuj k \R\ 2 
smh 2 (uj k T /2) 



[cosh(wfeT/2) cos/cAcc — 1] 



This implicitly determines T in terms of energy E. 
The exponent W at the saddle point equals 



W 



dk\R\ 



cos fcAcc — 1 



sinh(w fe T/2) e ^ T / 2 + l 



(64) 



(65) 



The integral of the second term is infrared-divergent. 
Using the longitudinal size L as an infrared cutoff, 
we find that to logarithmic accuracy the integral is 
(2ir I ' g) ln(L/T). In the exponent of eq. (|52|l for the prob- 
ability, W is combined with twice the instanton action 
l)55|l. As a result, the dependence on L disappears. 

By a direct calculation, one can verify that the saddle- 
point expressions l|64|l and (|65|l are related: 



dW 



= -E 



(66) 



This relation is convenient if we want to restore W (up 
to a constant) from an already calculated E. 

For real Ace, the first integral in l|65|) rapidly converges 
in the ultraviolet. If we use the simple acoustic dispersion 
law H46J1 , this integral can be computed explicitly, and we 
obtain, to logarithmic accuracy, 



2tt it Ax 

In cosh — — 

9 T 



2tt L 
— In — 
9 * 



(67) 



This is the same expression as obtained in ref. [4| by com- 
puting the action of a certain periodic field. For an indi- 
vidual real Acc, l|67|) is indeed an adequate approximation 
to W, but we still need to integrate over Ax. This will 
be done by steepest descent, and we will see that on the 
corresponding (complex) saddle point the integral does 
not converge as rapidly. As a result, the saddle point 
cannot be thoroughly explored in the acoustic approxi- 
mation: we will need the more accurate dispersion law 
I®. 

Using eq. H65|) . we obtain the following saddle-point 
equation for Acc: 



iP - 



sin 



fcA 



dkk\R\ 2 , 
«, smh(w fc T/2) 



(68) 



which shows that the saddle-point Acc lies on the upper 
imaginary axis. Eq. 168fl can be cast into a form similar 
to eqs. (£>EJ>, H3>: 



P 
~2 



dkka k a k 



dkkb k b k 



(69) 



The integrals here are the total momenta of quasipar- 
ticles in the initial and final states. Thus, on the one 
hand, <JHH| 

expresses momentum conservation and, on the 
other, shows that tunneling occurs between states with 
momenta =pP/2. 

We are interested in the range of energies, for which 
the periods satisfy 



1 

gn 



2£ 



(70) 



(the second relation applies due to our choice of units 
with c s = 1). In this case, the saddle-point Acc is of the 
form 



Ax = -(T + S) 



(71) 
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with < 5 <C T. Using the corrected dispersion law 
(|48|l in the exponents and the acoustic lo^ = \k\ in the 
coefficient R, we take eq. Ij68(l to the form 



- = -/ rffcexp(--fc 3 C 2 T+-fc ( 5). 
2 g J 4 2 



(72) 



We see that without the cubic term in the exponent the 
integrand would not have the correct large-A: behavior. 
The integrand in (|72Jl has a maximum at k = k*, 



K = 



25 

3^ 2 r 



(73) 



Assuming that the maximum is sufficiently sharp (this 
can be confirmed a posteriori), and approximating the 
integrand near it with a Gaussian, we obtain 



m- + 0(mln-) 



(74) 



Here we have used P = 27m = Tr/g£,- The sharpness of 
the maximum and therefore the accuracy of this calcula- 
tion is controlled by the large ln(T/£). 

Now, comparing the expressions l|58|l for the energy 
and H69|) for the momentum, we see that the main differ- 
ence between the two is due to the deviation of Uk from 
the strict acoustic form. Indeed, (assuming — > 0) we 
can write 



E - 



P 



dk(u k - k)b* k b k 



(75) 



Using (|61fl with the saddle-point values of C,, At, and Ax, 
we see that for k < the integral is rapidly converging, 
and the contribution from this region is small: most of 
the total comes from k > 0, where the integral converges 
much more slowly. Using the same Gaussian approxima- 
tion as above, we obtain 



E ~i = \eklp\l+0{\n^%) 



or, substituting fc* from lj73jl. 



E 



P 
~2 



P 

7 



21 

3T 



2/3 



In 



2/3 



T 

T 



(76) 



(77) 



Thus, the energy of the optimal initial state is close to 
P/2, but there is a correction, given by the right-hand 
side of {29 ■ 

Note that in the above calculation it was sufficient to 
use R(k) obtained in the limit u>k = \k\: any corrections 
to R(k) due to the modified dispersion law multiply the 
already small u>k — k in l|75|l and do not affect the leading 
correction computed in (|7(j[) . The same applies to any 
changes in the relation between a, b and R that are due 
to scattering corrections. 

We can now use (166(1 to restore the exponent W gov- 
erning the QPS rate. In doing so, we need to take into 



account the fact that in (|66|l the derivative is at fixed 
Ax, while l)76|) was obtained using the saddle-point Ax, 
which itself is a function of T . This difficulty can be 
circumvented in the following way. We first rewrite l)66|) 
as 



d_ 
df 



{W + iPAx)\ Ax = -E 



(78) 



and then observe that the partial derivative here can be 
replaced by the total, since the part due to the depen- 
dence of Ace on T vanishes at the saddle point. Thus, 
integrating (|76|l over T we obtain not just W but the 
sum W + iPAx. 

For the full exponent 



W = -2^0 + iPAx -iE( + W 



(79) 



which according to eq. H52[) governs the exponential fac- 
tor in the rate (in the limit En — > 0), we find 



P. 



7T /ST 



1/3 



T 



(80) 

where In = ln(T/£). In the same approximation, we can 
use (|77|l to express the period T through energy: 



p 



8 \E- P/2 



3/2 



hi 



P 



E -P/2 



(81) 



Substituting this into eq. H8()(l . we obtain our final ex- 
pression for the exponent: 



W » - — 



P 



1/2 



In- 



P 



(82) 



4g \E~P/2j E - P/2 ' 
which applies at energies such that 

£-P/2«P. (83) 
The microcanonical rate is 



U ~ exp(W) 



(84) 



It corresponds to the optimal choice of the tunneling end- 
point among all states of a given energy E. It is exponen- 
tially larger that the estimate l|2t)|) , which corresponds to 
some a priori, non-optimal, way of injecting the energy. 
In particular, the threshold for quasiparticle production 
has moved from E = P in to E = P/2 in 



D. The initial and final states 

We have referred several times to T as the period, but 
we have not yet exhibited the periodicity of the field con- 
figuration. The field can be obtained by substituting the 
saddle-point values of a and b* into eq. I|49|) . In addition 
to determining the field, these parameters also determine 
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the most probable initial and final states. We now turn 
to a detailed discussion of a and b*. 
The requisite saddle-point values are 



1 



a u)k At— ikAx 



K = 



1 — e Wk( 



iuikC, _ „uikAT — ikAx 



-Re 



(86) 



The values for a* and b are obtained by reverting all signs 
in all exponentials and replacing R with R* . Substituting 
these expressions into eq. ()49|l and setting v' — 0, we 
obtain the fluctuation v. 

The saddle-point solution 1)63(1 fixes the difference be- 
tween ro and Tq, but not the "center-of-mass" position 
(tq + To)/2. The latter will in general have both real and 
imaginary parts. The real part is fixed by the position 
of the time contour, on which the 5-matrix is defined. If 
we use the contour shown in Fig. [31 then the real part of 
To is zero, so we can write To = Hr, where tR is real — 
it is the moment (in real time) at which the QPS tales 
place. Since the collective mode associated with t R is 
not important for the present argument, we set tR = 0, 
so that 



(87) 



to = 0, 
r' Q = T/2 



This places the instanton at the origin, as indicated in 
Fig. El 

With the help of (|8"5)l and (|EEJl, the fluctuation v is 
obtained as an integral over k, see eq. ((4911 . To similarly 
represent the full field 9 = 9i + ^fgv, we need also the 
Fourier transform of the (modified) instanton field 8j. 
Note that this has different forms in the regions t > and 
t < 0, cf. eqs. lf4*4"|) . For definiteness, we consider 

t > 0, where we can use eq. ((44(1 with Tf replaced by t. 
As a result, the field can be written as 



-=9{X,T): 

dk 



iu) k C, 



G cjfc At — ikAx 



)/(as,r)E(A) , (89) 



where 



f(x t) = J^Q^kir-T^-ikix-xa) _|_ p-^k {r-T^)+ik(x-x' ) 

(90) 

and S(fc) is defined in ((54(1 . Note that S can be rewritten 
as 



E(fc)=£. 



(91) 



n=0 



This allows us to interpret ((89(1 as a sum over instantons 
and antiinstantons at various locations, in parallel with 
the interpretation of periodic instantons in ref. |19| . 

Indeed, if we substitute Q = iT and restrict our atten- 
tion to the interval (0, T/2), we see that the field 1(89(1 



can be interpreted as the sum of fields from two periodic 
chains: one of instantons, at locations x = xq, t = nT, 
and the other of anti-instantons, at locations x = x' , 
t = (?i + l/2)T, n = 0,±l,.... 

If we keep xo and x' real, we have a family of periodic 
configurations, such as the one shown in Fig. [I] These 
are the same configurations as found in ref. 4| by im- 
posing from the start the requirement of periodicity in 
the Euclidean time t. Here we have reconstructed them 
without any such a priori requirement, following instead 
the perturbative method of ref. We have seen, how- 
ever, that the integral over Ax in the probability l(52|) is 
not determine d by a real saddle point. So, unlike the case 
considered in [l!j , none of these real-zo configurations is 
an approximate classical solution. The approximate solu- 
tion that determines the rate corresponds to the complex 
saddle-point l(71|) and is itself complex. 

For real Xq, x' , the integral in l(89|) converges in the 
ultraviolet, for any r in the interval (0, T/2), even in the 
acoustic approximation lu^ = \k\. Using the acoustic dis- 
persion is equivalent to approximating each instanton in 
the chain by 9j(x — xq, t — t„), where 9i is the unper- 
turbed solution ((39(1 . The sum over n can then be done 
explicitly, resulting in |4( 



1 



( Xl T) = -[Hl-e- lT - x+x ")-Hl- 

1% 



ir — x-\-XQ \ 



- ln(l + e -**-*+*o) + ln(l + e iT - x+x °)] , (92) 

where all distances and times are measured in units of 
T/27T. This configuration is explicitly periodic with pe- 
riod 27r, and its Euclidean action per period to logarith- 
mic accuracy equals —W ac — iPAx + 2Sq, where W ac is 
the approximate expression l(67|) . 
Now, consider the case when 



XQ 



XR 
XR 



Ax/2 
Ax/2 



(93) 
(94) 



where xr is real, and Aa; is the saddle-point value 1(71(1 . 
Note that these Xq and x' are complex conjugate to each 
other. In this case, the integral in 1(89(1 is not convergent 
in the acoustic approximation, and we need to use the 
more precise dispersion law 1(48(1 . 

The sections of Euclidean time at t — ±T/4, which 
are half-way between the instantons and antiinstantons, 
are expected to have a special significance. Since the 
complete periodic solution determines the rate (i.e., the 
amplitude squared), we should be able to cut it in half, 
at t = ±T/4, to obtain the tunneling path, and then 
attach this tunneling path to real-time evolution. A di- 
rect calculation shows that when xq and x' are complex 
conjugate, as in l(9^1) . (|Mjl . the Euclidean velocity d T 9 at 
t = T/4 is purely imaginary, while the field itself there 
is purely real. The same is true at r — —T/4. Therefore, 
the solution becomes purely real at 



t = ±T/4 + it' 



(95) 
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with real t', i.e., on the horizontal, real-time, segments of 
the contour of Fig. [3J 

An immediate consequence of this result is that the 
real-time segments do not contribute to the imaginary 
part of the action, so the tunneling exponent is deter- 
mined by the Euclidean segment alone. Another con- 
sequence is that the real-time evolution can be inter- 
preted as formation and decay of the coherent fields cor- 
responding to the tunneling endpoints. For example, at 
r = T/4 + it' the real-time solution, as a function of t' 
and x, can be read off the expression (|89|) . in which we 
substitute eqs. I|93l) . I|94l) and the saddle-point values of 
all the parameters. Note that due to the deviation of the 
dispersion law from the purely acoustic one, the wave 
packets of quasiparticles do disperse, so at t' — ► ±00 the 
coherent states become collections of free quasiparticles. 
This means that the expressions (|58fl for the total energy 
and (|69[) for the total momentum should apply quite gen- 
erally, i.e., even with the scattering corrections included, 
despite the fact that these corrections will modify the 
simple relations l|85|) . I|86|) between a, b* and the instan- 
ton's Fourier transform. 

Finally, we note that, as clear from the expressions 
(|60|l . (|61|l for the quasiparticle densities, the initial state 
(at t' — > —00) contains mostly phonons with positive 
momenta, while the final state (at t' — > 00) mostly with 
negative. 



E. Canonical rate 

The leading exponential factor in the QPS rate at a 
temperature T = 1/(3 is obtained by integrating the 
microcanonical rate over energy E with the Boltzmann 
weight: 



where 



dEe-^TZ 



Wp = W-(3E 



dEe w ? 



(96) 



(97) 



Since all variational parameters in the expression 179|) 
for W are already at their saddle-point values, the to- 
tal derivative of W with respect to E coincides with the 
partial, i.e., 



dW 
~d~E 



T 



(98) 



Thus, the exponential in l|96(l quite generically has an 
extremum at the energy E* corresponding to T = (3. 
This is consistent with the standard argument, accord- 
ing to which at a finite temperature we should be look- 
ing for solutions that are periodic in r with period (3. 
However, in general, the extremum can be either a max- 
imum or a minimum. In the present case, under the 
low-temperature condition (3 3> £, the energy E* falls 



into the range (|83|l . where we can apply eq. (|82|l . We 
find that 



E* » - + - 



2/3 



In 2 / 3 | 



(99) 



is a maximum, and the only one in this range. Thus, 
levels with energies near E* give the main contribution to 
thermally-assisted tunneling among all levels in the range 
(|83|l . This has to compete with transitions (tunneling or 
over-barrier) from levels outside the range (JS3J), i.e., those 
for which E — P/2 = 0(P). However, the transition rate 
for such levels is suppressed at least by the Boltzmann 
factor exp(-E/T) = exp(-cP/T), where c > 1/2. On 
the other hand, as we will soon see, tunneling from E ps 
E„ is suppressed, at (3 3> £, mainly by exp(— P/2T), so 
at these temperatures it is more important. 

All that remains, then, is to substitute T — (3 = 1/T 
into eq. lj5U|l . We obtain 



(100) 



exp ■ 



^-^(^) 1/3 [l + 0(ln 



(3P 
Y 



3/3 
2? 



,-1/2 




which is our final result. Here In = ln(/3/£); the preex- 
ponent is estimated in the next subsection. 

Recall that P — 27m, where n is the average density. 
So, under the condition (|27|l . the first term in the expo- 
nent of (|10Ufl is much larger than the second. This term 
was obtained in ref. |J| from the approximate expression 
(|67|l . which is based upon using the acoustic dispersion 
law throughout. Obtaining the second, correction, term 
requires, as we have seen, the use of the more precise dis- 
persion law l|48(l . The second term becomes of the same 
order as the first at (3 ~ £, where the present approxima- 
tion breaks down. 



F. Estimate of the preexponent 

A naive dimensional estimate for the preexponent in 
the QPS rate would be L/£ 2 , where L is the total length 
of the system. In actuality, of course, there is a multi- 
plicative correction to this estimate. Note, however, that 
the estimate of accuracy in (|100fl implies that we have 
already omitted a variety of subleading terms, such as, 
for instance, 

2tt L 2tt (3 . 
- 2S Q + — In - « In £ 101 

9 P 9 ? 

(which is the leading contribution in the case of a sharply 
localized perturbation, but is only subleading in the uni- 
form system). At small g, this term, omitted in the expo- 
nent of the rate, is more important than any correction 
to the naive preexponent that we may obtain. Never- 
theless, we now present an estimate of the preexponent, 
given the traditional interest in values of the "attempt 
frequency" (which is what the preexponent represents). 
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The preexponent comes from the ratio of two deter- 
minants: one for small fluctuations near the periodic in- 
stanton, and the other for small fluctuations in vacuum — 
similarly to the case of a bounce |2^|, for a review, see ref. 
|28j. Near the periodic instanton, fluctuations can be ex- 
panded into normal modes, which satisfy the eigenvalue 
equation 



(102) 



where Dp is a small-fluctuation operator in the periodic- 
instanton background. The modes ip n should be peri- 
odic with the same period (3. We can use a classification 
of modes similar to that used after eq. I|51ll : there are 
modes that are localized in regions of a spatial size of 
order £ around the individual instantons and antiinstan- 
tons, and modes that are not. Delocalized modes with 
Ki <C l/^ 2 (i-e., with typical momenta much smaller than 
l/£), as well as modes with A„ > l/£ 2 , both delocalized 
and localized (if any), contribute only a numerical factor 
of order unity to the ratio of the determinants (cf. ref. 
|28jL and we will not be interested in this factor in what 
follows. 

Large factors can only come from localized modes with 
A,i <C l/£ 2 , corresponding to the "soft" collective coordi- 
nates. The periodic instanton has two strictly zero modes 
(A„ = 0), corresponding to translations of the entire con- 
figuration in space and time. These are the translations 
described by the parameters xr and tn of subsect. IIVDI 
Each of these zero modes contributes a normalization fac- 
tor of order up to a power of ln(/3/£), times the 
total volume associated with the corresponding collective 
coordinate. The temporal volume cancels out when we go 
from the probability to the rate, while the spatial volume 
L remains, resulting in the following zero-mode factor 



L 

W 2 



(103) 



In addition, there are two quasi-zero modes, correspond- 
ing to changes in the relative position of the instanton and 
antiinstanton chains. These modes are described by the 
saddle-point parameters Aa; and At. They have approx- 
imately the same normalization factors as the strictly 
zero modes, but their volumes are determined by the 
saddle-point integrations, i.e., by the second derivatives 
of the exponent W, eq. J53J), with respect to Ax and 
At. These volumes are of order (fc*P) -1 / 2 each (while 
the mixed second derivative vanishes), so the quasi-zero 
modes contribute the factor 



« 2 ~ (gfhP)- 1 ~ (tih)- 1 



(104) 



where k„ is given by eq. I|74|) with T = /3. Multiplying 
f2i and 0,2, we obtain an estimate for the preexponent: 



9 e u 



1/3 



(105) 



V. CONCLUSION 

In a uniform system, instantons that generate momen- 
tum by unwinding a persistent current need to transfer 
a compensating momentum to quasiparticles. We have 
seen that these instantons have many interesting prop- 
erties that are absent in cases when no such momentum 
transfer is necessary. We have considered in detail the 
case of a weakly-coupled ID superfluid at temperatures 
T -C gn ~ c s /£, where c s is the phonon speed, and £ is 
the healing length. (In this section, we restore c s , which 
was set to 1 earlier in the paper.), 

On the theoretical side, perhaps the most curious fea- 
tures of this case are that the instanton is complex, has no 
turning points, and yet its analytical continuation to the 
appropriate real-time segments is real. On the real-time 
segments, the solution can be interpreted as formation 
and decay of coherent states of quasiparticles. These ini- 
tial and final states have opposite total momenta =pP/2, 
allowing for a transfer of momentum P. 

A possibility of experimental detection of QPS in nar- 
row superfluid channels via momentum imaging has been 
mentioned in the introduction. Looking at our final re- 
sult (|100|) , which to the leading order we can rewrite as 



. 7rc s gn 
TZt ~ exp — 



g 



T 



(106) 



and comparing it to its counterpart p.'ill for the case of 
a sharply localized perturbation, we see that the tem- 
perature dependence of the rate is much steeper in the 
uniform case: a logarithm of gn/T in the exponent is now 
replaced by a power. As the calculation makes clear, this 
additional suppression results from the need to use states 
with relatively high energy E ss c s P/2. The "dynamical" 
suppression, controlled by the quantum overlap between 
the optimal initial and final states, is in the uniform case 
only subleading. 

It is of interest to extend the present results in several 
directions. First, it would be interesting to extend them 
to higher temperatures, T ~ gn. A striking property 
of the rate (|106|) is that, at any T <C gn, it is exponen- 
tially smaller than the rate of thermal activation that one 
would obtain within the LAMH theory [l(J, [Tl| . Indeed, 
when the energy Ejq released by unwinding the current is 
negligible (the same limit, in which (|106|) was obtained), 
the LAMH rate is 



"ft-LAMH ~ exp 



4c s n 
3T 



(107) 



which is accurate up to a power of ln(/3/£). 



On the other hand, we have shown in subsect. IIV El that 
at T <C gn our instanton is the dominant path for phase 
slips, more important than any thermal activation. The 
discrepancy can be traced to the fact that the original 
LAMH theory makes no account of momentum conser- 
vation. Recall that the LAMH saddle point, whose en- 
ergy determines the rate (|107H , has order parameter that 
is nonvanishing everywhere (except for the special case 
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of exactly zero current). Usually, one assumes that this 
saddle point is close to some time-dependent fluctuation, 
for which the order parameter vanishes at some point, 
allowing for a phase slip. Our results imply that in a 
uniform system, where there are no external "sinks" of 
momentum (such as impurities, etc.), this is not a good 
assumption, i.e., in this case the LAMH saddle point does 
not nucleate any phase-slip process. We have preliminary 
numerical data supporting this conclusion. Accordingly, 
it is by no means clear if, in a uniform system, there is a 
crossover to a thermally activated mechanism for phase 
slips at any T . 

Another direction in which one may be able to extend 
the present results is systems of higher dimensionality 
(2D and 3D). There, the role of the topological term is 
taken over by the Magnus force, which suppresses tun- 
neling in a rather similar way [LJ]. In either case, the 
suppression can be seen as a result of destructive interfer- 
ence between tunneling at different values of the spatial 



coordinate. It is natural to ask if in 2D and 3D this sup- 
pression can be circumvented by an inelastic mechanism 
(production of phonons) similar to the one considered 
here. 

One last case of interest we mention is that of BCS- 
paired superfluids and superconductors. In that case, an 
additional channel of momentum production £| becomes 
available. It is associated with zero modes of fermionic 
quasiparticles at the instanton core. For superconduc- 
tors, the problem is complicated by the scattering of 
quasiparticles on disorder (and, in ID and 2D, on the 
boundaries of the sample), which alters the momentum 
balance. Nevertheless, this channel of inelastic tunneling 
deserves a further study, especially in view of its possi- 
ble relevance to experiments |lo . Il7j on superconducting 
nanowires. 
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